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Abstract 

A summary of a new approach is given for analysing the 
dynamic behaviour of distributions of charged particles in 
an electromagnetic field. Noting the limitations inherent 
in the Lorentz-Dirac equation for a single point particle, a 
simple model is proposed for a charged continuum inter¬ 
acting self-consistently with the Maxwell field in vacuo. 
The model is developed using intrinsic tensor field the¬ 
ory and exploits to the full the symmetry and light-cone 
structure of Minkowski spacetime. A covariant perturba¬ 
tion scheme is motivated by an exact class of solutions to 
this system describing the evolution of a charged fluid un¬ 
der the combined effects of both self and external electro¬ 
magnetic fields. The scheme yields an asymptotic approx¬ 
imation involving inhomogeneous linear equations for the 
self-consistent Maxwell field, charge current and time-like 
velocity field of the charged fluid 

INTRODUCTION 

The intense international activity involved in probing the 
structure of matter on all scales, with particle beams and 
radiation, owes much to recent advances in accelerator sci¬ 
ence and technology. Developments in the production of 
high power laser radiation also offer new avenues for accel¬ 
erator design and new diagnostic tools of relevance to med¬ 
ical science, engineering and the communications industry. 
A common theme in these developments is the interaction 
between charged particles and the electromagnetic field in 
domains where relativistic effects cannot be ignored. 

It is remarkable that many of the challenges that must 
be addressed in order to develop and control devices that 
accelerate charged particles have their origin in the interac¬ 
tion of particles with their own electromagnetic field. De¬ 
spite the fact that the classical laws of electromagnetism 
were essentially formulated over a century and a half ago 
the subject of electromagnetic interactions with matter re¬ 
mains incomplete. This incompleteness has had concomi¬ 
tant effects on the development of quantum electrodynam¬ 
ics and renormalisation theory. At root, the difficulties re¬ 
side in the recognition that the quantum structure of matter 
at some scale is beyond observation. Furthermore, the clas¬ 
sical description of the electron as a point particle leads to 
singularities in the Maxwell self-fields that inevitably cre¬ 
ate ambiguities in its interaction with the Maxwell field. 
The general consensus is that a useful domain of validity of 
the Lorentz-Dirac equation [1,2, 3], describing covariantly 
the radiation reaction on a point electron, can be accommo¬ 
dated by performing a “reduction of order” that effectively 
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replaces the equation by a perturbative second order system 
for the particle world line. One must then decide whether 
higher order terms in this expansion should be maintained 
given the neglect of terms associated with the regularisation 
scheme. This approach is behind many of the successful 
applications of approximate radiation reaction dynamics, 
despite the somewhat delicate and unsatisfactory nature of 
the arguments that purport to support it. However, an un¬ 
derstanding of the interaction of charged particle bunches 
with their electromagnetic fields is vital for the success of 
proposed future accelerators. 

In this article a summary is given of a new approach for 
analysing the behaviour of distributions of ultra-relativistic 
charged particles in a coupled electromagnetic field en¬ 
vironment [4]. Rather than mixing particle and field 
concepts, our model is formulated entirely in the lan¬ 
guage of smooth tensor fields. Here, the non-linear field 
equations describing a thermodynamically inert (pressure¬ 
less) charged fluid are adopted as a dynamical model of 
charged particle bunches. They are analysed using an ultra- 
relativistic asymptotic approximation scheme that con¬ 
serves total energy, momentum and charge at all orders in 
the expansion parameter. 


EQUATIONS OF MOTION FOR A 
PRESSURELESS CHARGED FLUID 

The partial differential system of equations governing 
the motion of a pressureless charged fluid in MKS form 
is 

V • e = j-g, V x b = /j, 0 qv + §f, 

Vxe+§=0, V • 6 = 0, 
fE + ( v .V)p=«, (e + uxb), 



where e is the electric field, b is the magnetic induction, v 
is the fluid’s velocity field, g is the fluid’s charge density, 
p is the fluid’s momentum field, mo is the rest mass of the 
particles comprising the fluid, q 0 is the fundamental charge 
of the particles comprising the fluid, eo is the permittivity 
of the vacuum and po is the permeability of the vacuum. 
Although (1) is the most familiar statement of the fluid’s 
equations of motion, their most succinct and analytically 
powerful form is expressed without reference to a particular 
inertial frame i.e. where the inertial time t does not appear 
explicitly. Let M be Minkowski spacetime and let g be its 
flat Lorentzian metric with V the Levi-Civita connection. 
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The system (1) may be recast in the coordinate-free form 

dF = 0, d*F = -p*V, 

V v V = i v F, g(V,V) = -1 

where V is the fluid’s 4-velocity field, F is an electromag¬ 
netic field 2-form, p is proportional to the proper charge 
density 1 and is called the reduced proper charge density, 
d is the exterior derivative, * is the Hodge map associated 
with g, iy is the interior operator and the 1-form V is the 
metric dual of V. For more information on the precise re¬ 
lationship between (2) and (1) see [4], 

APPROXIMATION SCHEME 

In general situations, exact solutions to the non-linear 
partial differential system (2) are extremely difficult to find. 
A self-consistent method for obtaining approximate solu¬ 
tions is to let V, F and p depend on a running parameter e 
and to solve the hierarchy of equations generated by (2) for 
the 1-parameter families of fields V s , F e and p e order by 
order in e. Whether or not the hierarchy of equations so ob¬ 
tained is easier to solve than the original system depends on 
the £ dependences of the fields. Although there are many 
possibilities, a viable choice is strongly suggested by the 
behaviour of a class of highly symmetric exact solutions 
relevant to accelerator physics. These solutions describe 
“walls of charge” of infinite extent accelerating due to their 
self-field (“space-charge”) and an applied constant electric 
field. Let (t, x, y, z) be an inertial coordinate system where 
t is the laboratory time and (x. y. z) are Cartesian coordi¬ 
nates. The spacetime metric tensor g has the form 2 

g = —dt <g) dt + dx® dx + dy ® dy + dz ® dz. (3) 

The “wall of charge” solutions have the following proper¬ 
ties in the (dt, dx, dy, dz) frame : the magnetic field van¬ 
ishes, the electric field has only a 2 component, the 4- 
velocity field has only t and z components and all field 
components depend only on t and z. The z component 
of the reduced electric field has the form 

£(t,z) = C,(d(t,z)) (4) 

where a = d(t, z) is a solution to the implicit equation 

( 5 ) 

and the lines of constant a are the integral curves of V. 
Splitting the initial (f = 0) electric field into a constant 
external part C,~\/e and a self part (o(z)\ 

£ s (0,z) = ^ t + <;o(z) ( 6 ) 

'The proper charge density is e °™° c P■ 

2 Units are chosen in which the speed of light c = 1. 


leads to 

£ £ (t, z) = jC-i + Co(2 — t) + + 0(s 2 ), 

for f > 0 and C-i > 0 and where C 0 ((t) = ^(<7). Using 
(7), it can be shown that the 4-velocity field V e and reduced 
proper charge density p e take the form 

V £ = [jC-l + Co(^ — f)j t (d t + d z ) 

•fr p +2 ggf dt + s + 0(e 2 ), ( g ) 

p £ = s^ + 0(e 2 ). 

Note that the coefficient of 1/e in V s is light-like; there¬ 
fore, such field configurations are called ultra-relativistic. 
More generally, it can be shown that the £ expansions of 
F s , V s and p e take the form 

F e = £ e n F n , V e 

P S = 

where F_ 1 is an external field (a solution to the source-free 
Maxwell equations). Therefore, (9) is adopted as an ansatz 
for the £ dependence of general solutions that include the 
effects of guiding and accelerating electromagnetic fields. 
Indeed, inserting (9) in (2) leads to a partially decoupled (in 
comparison with (2)) infinite set of differential equations 
whose myriad solutions include (7) and (8). 

The infinite system is solved order by order in e accord¬ 
ing to the following iterated pattern: 

solve Maxwell’s equations 

I 

solve for the motion of the continuum (10) 

I 

enforce charge conservation 

The first step in the solution procedure is to note that the ex¬ 
ternal electromagnetic field F_i is a solution to the source- 
free Maxwell equations: 

dF_ i=0, d* F_i = 0. (11) 

The second step involves solving 

Vv- 1 V- l = i v _ 1 F_ 1 , g(V- 1 ,V- 1 ) = 0 (12) 

for the leading order 4-velocity field U_-|. Note that, to 

leading order, the motion of the charged continuum is gov¬ 
erned by the external field F_ 1. The third step is to enforce 
charge conservation to leading order by solving 

d*(piU_i)=0 (13) 
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for pi and ensures that the following step is consistent. The 
fourth step is to solve the Maxwell equations for Fq with 
the 4-current piV_i as a source: 

dFo = 0, d*F 0 = -k Pl V- 1 (14) 


and the fifth step is to solve 


Vy_i Vo + Vy 0 E_i = i V _ 1 Fo + iy 0 F_i, 

g(V-i,V 0 )= 0 


(15) 


for Vo. The sixth step is to enforce charge conservation to 
next-to-leading order by solving 

d* (p 2 V-i) + d* (piE 0 ) = 0 (16) 

for p 2 , leading to the seventh step which is to solve the 
Maxwell equations 

dF 1 = 0, d*Fi = — *p 2 V_ 1 - * Pl V 0 (17) 


for Ei. This procedure is continued to any order in e and 
leads to expressions for F e , V s and p e to any desired level 
of accuracy. 

The key point to note is that the only equation non-linear 
in its unknown (the vector field V-i) is (12). However, (12) 
is straightforward to analyse because it can be reduced to a 
quasi-linear second order ordinary differential equation for 
the integral curves of IF j . It can be shown that the equa¬ 
tions for V n with n > 0, such as (15) for Vo, written out 
in a frame adapted to V_i lead to to three inhomogeneous 
linear differential equations determining three of the com¬ 
ponents of V n . The remaining component is the solution to 
a linear algebraic equation involving fields of lower order 
in s. 


EXAMPLE 


where #±1 and dj_ are the volume 2-form and exterior 
derivative in the ( x , y) planes. The 3-velocity of the beam 
is along the z axis and has the subluminal Newtonian speed 

1 - ^ + 0(£ 3 )- 

A Gaussian bunch with transverse radius Ro travelling 
at constant Newtonian speed 1 — + 0(s 2 ) is described 

by 

7-i {x,y) = b 0 , 

Pi(z,x,y) = a 0 exp (~^-) S(*) (21) 

where ao, Ro and l>o are constants and S is a smooth 
bump function vanishing outside the interval (—z±,zi) and 
S(z) = 1 for z € (-Z2, z 2 ) and z± > z 2 > 0. The leading 
order laboratory reduced charge density 7-ipi for some 
range of t is 


7-i Pi =l-i(x,y)pi(z-t,x,y) 

= a 0 b 0 exp x +■? ) E(z-t). ^ 

Working in the cylindrical polar coordinates (t, R, <p, z) 
where x = R cos <j> and y = R sin (p, a cylindrically sym¬ 
metric solution to (20) well-behaved at R = 0 is 


4, ° = {i [ 1 -“p(-^)]*}^-<) 

(23) 

and the corresponding electromagnetic 2-form F t) is 

F ° = a ° b °W, f 1 “ exp ( _ ^)] dRA (~ dt + dz )- 

(24) 

The laboratory electric field is radial, the magnetic field is 
azimuthal and their magnitudes are equal and vanish out¬ 
side of the support of H. 


A simple application of the approximation scheme is to 
consider a high-energy charged beam propagating in free 
space where the external field F_i is zero. Following the 
solution method introduced in the previous section, it may 
be shown that the following 1-parameter fields are approx¬ 
imate solutions to (2): 

^ £ = (b-1 + i 

%||7-i 4 d z + 0(e 2 ) (18) 

F e = — (d $ 0 + ed$i) A dt 

+(d $o + ed<f>i) Adz + 0 (s 2 ) 

where 

7-1 = 7 (x,y), 

Pi = Pi{z - t, x, y), p 2 = p 2 (z - t, x, y) 
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are fixed by by their forms at an instant in t given as data. 
The fields 4>o and 4> i are solutions to the transverse Poisson 
equations 

d±#±d±$o = 7—iPi#-lV r9m 

d±#±d±$r=7- 1 P2#±l ( ’ 
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